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The carrier of the "30" /xm emission feature in evolved stars^ 

A simple model using magnesium sulfide 
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Abstract. We present 2—45 /jm spectra of a large sample of carbon-rich evolved stars in order to study the "30" 
(j,m feature. We find the "30" /im feature in sources in a wide range of sources: low mass loss carbon stars, 
extreme carbon-stars, post-AGB objects and planetary nebulae. We extract the profiles from the sources by using 
a simple systematic approach to model the continuum. We find large variations in the wavelength and width of 
the extracted profiles of the "30" /im feature. We modelled the whole range of profiles in a simple way by using 
magnesium sulfide (MgS) dust grains with a MgS grain temperature different from the continuum temperature. 
The systematic change in peak positions can be explained by cooling of MgS grains as the star evolves off the 
AGB. In several sources we find that a residual emission excess at ~26 fjxa can also be fitted using MgS grains but 
with a different grains shape distribution. The profiles of the "30" fim feature in planetary nebulae are narrower 
than our simple MgS model predicts. We discuss the possible reasons for this difference. We find a sample of warm 
carbon-stars with very cold MgS grains. We discuss possible causes for this phenomenon. We find no evidence 
for rapid destruction of MgS during the planetary nebula phase and conclude that the MgS may survive to be 
incorporated in the ISM. 
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>■ . 1. Introduction 



1985; Waters et al. 200C; Hrivnak et al. 2000; Volk et al 



The far infrared (iKj spectra ot carbon-rich evolved ob- 
jects; i.e., carbon-rich AGB stars (C-stars), post asymp- 
totic giant branch objects (post-AGBs) and planetary 
nebulae (PNe) are typified by a broad emission feature 
around 30 (ira. This "30" fim feature was first discovered 
in the far-IR spectra of CW Leo, IC 418 and NGC 6572 by 
Forrest et al. (1981). Since then this feature has been de 



2002) 



Goebel & Moseley (1985) proposed solid magnesium 
sulfide (MgS) as the possible carrier of the "30" /im fea- 
ture. Their suggestion is based on the coincidence of the 
emission feature with the sole IR-resonance of MgS (Nuth 
|et al. 1985| ; [Begemann et al. 1994] ) and the fact that MgS 
is one of the expected condensates around these objects 



tectedm" a wid7ranW"carbori-rich evolved objectTfrom (^timcr et aL 197 ^ fodders fc Fegley 19990. Several 



intermediate mass loss C-stars (Yamamura et al. 1998) 
to post-AGBs and PNe (pmont 1993j; pox 1993 



et al. 1995; Jiang et al. 1999; Bzczerba et al. 1999; Hony 



Omont 



authors have taken up on this suggestion and compared 
observations with laboratory measurements of MgS. These 
comparisons were further facilitated by the publication of 



et al. 2001 ). The feature is commonly found in C-rich post- 
AGBs and PNe however with varyin g band shapes and 
varying feature to continuum ratios ( Goebel fc Moseley 



the optical constants of MgS in the IR range by Begemann 
|ct al. (1994 ). These authors found that the far IR excess 
of CW Leo can be successfully modelled using MgS grains 
with a broad shape distribution. 
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based on observations obtained with ISO, an ESA project 



with instruments funded by ESA Member states (especially 



the PI countries: France, Germany, the Netherlands and the 



United Kingdom) with the participation of ISAS and NASA 



More recently, Jiang et al. (1999) and Szczerba et al 
(1999| ) have modelled the spectra taken with the Short 



Wavelength Spectrometer (SWS) ( de Graauw et al. 199G ) 
on-board the Infrared Space Observatory (ISO) (Kesslei 
|ct al. 1996Q of the C-star IRAS 03313+6058 and the post- 
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AGB object IRAS 04296+3429 respectively. They find 
that for these sources which show a strong "30" /im fea- 
ture, the elemental abundances of Mg and S are consistent 
with MgS as the carrier of the feature. 
' Hrivnak et al. (200fj) and |Volk et al. (2002[ ) have anal- 
ysed ISO spectra of a sample of post-AGBs. They find 
that the profile of the "30" /im feature varies between 
sources. Although these authors state that this decompo- 
sition is not unique, they find that their "30" /im feature 
is composed of two sub features: one feature peaking near 
26 /im and an other near 30 /im. Using these two compo- 
nents in varying relative amounts they are able to explain 
the range of features found in their sample. Based on the 
discovery of these sub features they consider the carrier(s) 
of the "30" /im feature to be unidentified. 

Other materials have also been proposed as carriers 
of the "30" /im feature. Duley (2000 ) suggest that the 
"30" /xm feature may be indicative of carbon-based linear 
molecules with specific side groups. Such molecules have 
strong absorption bands throughout the 15—30 /im range. 
Papoular (2000 ) discusses the possible contribution of car- 
bonaceous dust grains with oxygen in the structure. Some 
of these materials may show IR emission in the 20—30 
/im range. Since the optical properties of such grains are 
sensitive to the exact composition they might be able to 
explain the range of features found in the C-rich evolved 
stars. Recently, Grishko et al. (2001) have proposed hydro- 
genated amorphous carbon (HAC) as a possible carrier of 
the "30" /mi feature. 

The ISO mission has provided an excellent database of 
observations to study the properties of the "30" /im fea- 
ture in detail and test the suggested identifications system- 
atically. The wavelength coverage of the SWS instrument 
(2-45 /im) is sufficient to determine a reliable continuum. 
The sensitivity of the ISO spectrograph allows detection 
of relatively weak features. The resolving power of the in- 
strument (A/AA = 500-1500) makes it feasible to study 
possible substructure in the "30" /im feature. Thus these 
observations allow a study of the "30" /im feature in un- 
precedented detail in a large sample of sources. 

In this paper, we investigate the shape and strength of 
the "30 /im" in a wide range of objects from visual visible 
C-stars, extreme C-stars, post-AGBs to PNe in order to 
further test the MgS or other identifications and map sys- 
tematic differences between the feature in different classes 
of sources. 

Our paper is organised as follows. In Sect. ||, we de- 
scribe the sample and the data reduction. In Sect. [| we 
present the way in which we modelled the continuum in or- 
der to extract the feature properties. In Sect. ^, we present 
the full range of extracted profile shapes and peak posi- 
tions of the "30" /xm feature and we discuss the possible 
ways of interpreting the observed profiles. In Sect. [|, we 
develop a simple model using MgS for the "30" /im fea- 
ture. In Sect. ^, we present the model results and compare 
them to the astronomical spectra. In Sect. |^, we present a 
correlation study between several feature properties and 
stellar parameters. Finally, in Sect. |L we discuss the im- 
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Fig. 1. The IRAS two-colour diagram for the sources stud- 
ied in this sample. The triangles represent the C-stars, the 
squares are the post-AGB objects, the stars are the PNe 
and the diamonds are the C-stars without a "30" /im fea- 
ture detected. The dashed line represents the position of 
a blackbodies of different temperatures. The dotted line 
sketches the evolution of a C-star with a detached, ex- 
panding and cooling circumstcllar shell. 

plications of our model results and the consequences for 
the MgS identification. In particular, we discuss possible 
causes for the deviating profiles and the possibility that 
MgS produced in carbon-rich evolved stars will be incor- 
porated in the interstellar medium (ISM). 

2. Observations 

We present observations obtained with ISO of a sample 
of bright IR sources at different stages along the evolu- 
tionary track from C-star via post-AGB object to PN. 
The observations presented here consist of data obtained 
with the ISO/SWS using astronomical observing template 
06 and 01 at various speeds. These observing modes pro- 
duces observations from 2.3 to 45 /im with a resolving 
power(A/AA) ranging from 500 to 1500. The sample con- 
sists of all carbon-rich evolved objects in the ISO archive 
which exhibit a "30" /im feature stronger than 8 Jy peak 
intensity and have been observed over the full 2.3—45.2 
/im wavelength range of the ISO/SWS. This peak intensity 
and the typical noise level of SWS band 4 (29—45.2 /im) 
allows to extract a reasonably reliable feature strength and 
profile. The complete wavelength coverage is needed in or- 
der to provide a sufficient baseline to estimate the contin- 
uum. We have further completed the sample with all ob- 
served C-stars with an IRAS 25 /im flux over 13 Jy. These 
sources serve as a control group since we would expect to 
detect the "30" /im feature based on this brightness, the 
typical noise levels and the typical feature over continuum 



Hony et al.: The carrier of the "30" fim feature 
Table 1. Source list. Observational details of the sources in this study. 
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Object 


IRAS name 


Obs." 


a 


<5 


TDT" 


Sp./T 


Obj.Type 






Mode 


(J2000) 


(J2000) 




kK 




NGC 40 


00102+7214 


01(3) 


001301.10 


+72 31 19.09 


30003803 


WC 


PN 


IRAS 00210+6221 


00210+6213 


01(1) 


00 23 51.20 


+62 38 07.01 


40401901 




C-star 


IRAS 01005+7910 


01005+7910 


01(2) 


0104 45.70 


+79 26 47.00 


68600302 


OBc 


post-AGB 


HV Cas 


01080+5327 


01(1) 


01 1103.50 


+53 43 40.30 


62902503 




C-star 


RAFGL 190 


01144+6658 


01(2) 


01 1751.60 


+6713 53.90 


68800128 




C-star 


R Scl f 


01246-3248 












C-star 






01(2) 


01 26 58.10 


-32 32 34.91 


37801213 










01(2) 


01 26 58.05 


-32 32 34.19 


37801443 






IRAS Z02229+6208 


Z02229+6208 


01(1) 


02 2641.80 


+62 2122.00 


44804704 


GO 


post-AGB 


RAFGL 341 


02293+5748 


01(1) 


02 33 00.16 


+58 02 04.99 


80002450 




C-star 


IRC+50 096 


03229+4721 


01(2) 


03 26 29.80 


+473147.10 


81002351 




C-star 


IRAS 03313+6058 


03313+6058 


01(1) 


03 35 31.50 


+6108 51.00 


62301907 




C-star 


U Cam 


03374+6229 


01(2) 


03 4148.16 


+62 38 55.21 


64001445 




C-star 


RAFGL 618 


04395+3601 


01(3) 


04 42 53.30 


+36 06 52.99 


68800561 


B0 


PN 


W Ori 


05028+0106 


01(3) 


05 05 23.70 


+01 10 39.22 


85801604 




C-star 


IC 418 


05251-1244 


01(2) 


05 2728.31 


-12 4148.19 


82901301 


36 1 


PN 


V636 Mon 


06226-0905 


01(1) 


06 25 01.60 


-09 0716.00 


86706617 




C-star 


RAFGL 940 


06238+0904 


01(2) 


06 26 37.30 


+09 02 16.01 


87102602 




C-star 


IRAS 06582+1507 


06582+1507 


01(2) 


070108.40 


+15 03 40.00 


71002102 




C-star 


HD 56126 f 


07134+1005 










F5 


post-AGB 






06 


071610.20 


+09 59 48.01 


71802201 










06 


071610.30 


+09 59 48.01 


72201702 










01(3) 


071610.20 


+09 59 48.01 


72201901 






CW Leo 


09451+1330 


06 


09 4757.27 


+1316 42.82 


19900101 




C-star 


NGC 3918 


11478-5654 


01(1) 


115018.91 


-5710 51.10 


29900201 




PN 


RU Vir 


12447+0425 


01(2) 


12 4718.43 


+04 08 41.89 


24601053 




C-star 


IRAS 13416-6243 


13416-6243 


01(3) 


13 45 07.61 


-62 58 18.98 


62803904 




post-AGB 


II Lup 


15194-5115 


06 


15 23 04.91 


-5125 59.02 


29700401 




C-star 


V Crb 


15477+3943 


06 


154931.21 


+39 3417.80 


25502252 




C-star 


PN K 2-16 f 


16416-2758 










WC 


PN 






01(1) 


16 44 49.10 


-28 04 05.02 


29302010 










01(2) 


16 44 49.10 


-28 04 05.02 


67501241 






IRAS 16594-4656 


16594-4656 


01(1) 


17 03 09.67 


-4700 27.90 


45800441 




post-AGB 


NGC 6369 


17262-2343 


01(1) 


17 29 20.80 


-2345 32.00 


45601901 


WC8 2 


PN 


IRC+20 326 


17297+1747 


01(1) 


17 3154.90 


+1745 20.02 


81601210 




C-star 


CD-49 11554 


17311-4924 


01(2) 


17 35 02.41 


-49 26 22.31 


10300636 


BHIe 


post-AGB 


PN HB 5 


17447-2958 


01(3) 


174756.11 


-29 59 39.70 


49400104 




PN 


RAFGL 5416 


17534-3030 


01(1) 


17 56 36.90 


-30 30 47.02 


12102004 




C-star 


T Dra 


17556+5813 


01(2) 


17 56 23.30 


+5813 06.38 


34601702 




C-star 


RAFGL 2155 


18240+2326 


01(1) 


18 26 05.69 


+23 28 46.31 


47100261 




C-star 


IRAS 18240-0244 


18240-0244 


oi(i) 


18 26 40.00 


-02 42 56.99 


14900804 


WC 


PN 


IRC+00 365 


18398-0220 


01(2) 


18 42 24.68 


-02 1725.19 


49901342 




C-star 


RAFGL 2256 


18464-0656 


01(1) 


18 4910.35 


-06 53 03.41 


48300563 




C-star 


PN K 3-17 


18538+0703 


01(2) 


18 5618.05 


+0707 25.61 


49900640 




PN 


IRC+10 401 


19008+0726 


oi(i) 


19 0318.10 


+0730 43.99 


87201221 




C-star 


IRAS 19068+0544 


19068+0544 


01(1) 


19 0915.40 


+05 49 05.99 


47901374 




C-star 


NGC 6790 


19204+0124 


01(1) 


19 22 57.00 


+0130 46.51 


13401107 


70 3 


PN 


RAFGL 2392 


19248+0658 


01(1) 


19 2714.40 


+0704 09.98 


85800120 




C-star 


NGC 6826 


19434+5024 


01(4) 


19 44 48.20 


+50 3130.00 


27200786 


50 4 


PN 


IRAS 19454+2920 


19454+2920 


oi(i) 


19 4724.25 


+29 2811.78 


52601347 




post-AGB 


HD 187885 


19500-1709 


01(2) 


19 52 52.59 


-170149.58 


14400346 


F2 


post-AGB 


RAFGL 2477 


19548+3035 


01(1) 


19 56 48.26 


+30 43 59.20 


56100849 




C-star 


IRAS 19584+2652 


19584+2652 


01(1) 


20 00 31.00 


+2700 37.01 


52600868 




C-star 


IRAS 20000+3239 


20000+3239 


01(1) 


20 0159.50 


+32 47 33.00 


18500531 


G8 


post-AGB 


VCygt 


20396+4757 












C-star 






01(2) 


20 41 18.20 


+48 08 29.00 


42100111 










01(2) 


20 41 18.20 


+48 08 29.00 


42300307 
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Table 1. (continued). 
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Object 


IRAS name 


Obs. a 


a 


5 


TDT b 


Sp./T 


Obi. Type 






Mode 


(J2000) 


(J2000) 




kK 




NGC 7027 T 












200 b 


PN 






01(4) 


21 0701.71 


+42 14 09.10 


02401183 










01(1) 


21 0701.70 


+42 14 09.10 


23001356 










01(2) 


21 0701.70 


+42 14 09.10 


23001357 










01(3) 


21 0701.70 


+42 14 09.10 


23001358 










01(4) 


21 0701.63 


+42 14 10.28 


55800537 






S Cep 


21358+7823 


01(1) 


21 35 12.80 


+78 37 28.20 


56200926 




C-star 


RAFGL 2688 




01(3) 


21 02 18.80 


+36 41 37.79 


35102563 


F5 


post-AGB 


RAFGL 2699 


21027+5309 


01(1) 


21 04 14.70 


+53 21 02.99 


77800722 




C-star 


IC 5117 


21306+4422 


01(1) 


21 32 30.83 


+44 35 47.29 


36701824 


77 3 


PN 


RAFGL 5625 


21318+5631 


01(1) 


21 33 22.30 


+5644 39.80 


11101103 




C-star 


IRAS 21489+5301 


21489+5301 


01(1) 


21 50 45.00 


+5315 28.01 


15901205 




C-star 


SAO 34504 


22272+5435 


01(2) 


22 29 10.31 


+54 51 07.20 


26302115 


G5 


post-AGB 


IRAS 22303+5950 


22303+5950 


01(1) 


22 32 12.80 


+60 06 04.00 


77900836 




C-star 


IRAS 22574+6609 


22574+6609 


01(2) 


22 59 18.30 


+66 25 49.01 


39601910 




post-AGB 


RAFGL 3068 


23166+1655 


01(2) 


231912.48 


+171133.40 


37900867 




C-star 


RAFGL 3099 


23257+1038 


01(1) 


23 2816.90 


+10 54 40.00 


78200523 




C-star 


IRAS 23304+6147 


23304+6147 


01(3) 


23 32 44.94 


+62 03 49.61 


39601867 


G2 


post-AGB 


IRAS 23321+6545 


23321+6545 


01(1) 


23 34 22.53 


+66 0150.41 


25500248 




post-AGB 


IRC+40 540 


23320+4316 


01(2) 


23 34 27.86 


+43 33 00.40 


38201557 




C-star 


non detections 


R For 


02270-2619 


01(1) 


02 2915.30 


-26 05 56.18 


82001817 




C-star 


SS Vir 


12226+0102 


01(1) 


12 25 14.40 


+0046 10.20 


21100138 




C-star 


Y CVn 


12427+4542 


01(2) 


12 45 07.80 


+45 26 24.90 


16000926 




C-star 


RY Dra 


12544+6615 


01(3) 


12 56 25.70 


+65 59 39.01 


54300203 




C-star 


C* 2178 


14371-6233 


01(1) 


14 4102.50 


-62 45 54.00 


43600471 




C-star 


V1079 Sco 


17172-4020 


01(1) 


17 20 46.20 


-40 23 18.10 


46200776 




C-star 


T Lyr 


18306+3657 


06 


18 32 19.99 


+36 59 55.50 


36100832 




C-star 


S Set 


18476-0758 


01(2) 


18 5019.93 


-0754 26.39 


16401849 




C-star 


V Aql 


19017-0545 


01(2) 


19 04 24.07 


-05 4105.71 


16402151 




C-star 


V460 Cyg 


21399+3516 


01(1) 


214201.10 


+35 30 36.00 


74500512 




C-star 


PQ Cep 


21440+7324 


01(1) 


2144 28.80 


+73 38 03.01 


42602373 




C-star 


TX Psc 


23438+0312 


06 


23 46 23.57 


+03 29 13.70 


37501937 




C-star 



a SWS observing mode used (see dc Graauw ct al. 1996| ). Numbers in brackets correspond to the scanning Ipeed. 
^TDT number which uniquely identifies each ISO observation. 

' These spectra have been obt ained by co-adding t he separate SWS sp ec tra also listed in the ta bl e, see text. 
Effec tive temperatures f rom 1 Mendez et al. (1992), 
and 5 [Latter et al. (2000|). 



^Perinotto (1991), J Kaler fc Jacoby (1991), 4 Quigley fc Bruhweiler (1995) 



level. These sources without the "30" p feature detected 
are listed separately in Table |l| It should be emphasised 
that the ISO archive does not contain a statistically repre- 
sentative sample of objects. The database of observations 
for the carbon stars provides a reasonable sampling over 
stellar properties (e.g. mass-loss rates or colour temper- 
atures). However the post-AGB sample is heavily biased 
towards the "21" /im objects; a peculiar type of C-rich 
post-AGB object. The sample of PNe contains a collection 
of either bright, well-known or well-studied objects with- 
out a proper statistical selection. It also contains a rela- 
tively large proportion of PNe with hydrogen-poor central 
stars. The total sample of 75 sources contains 48 C-stars, 
14 post-AGB objects and 13 PNe. We have detected the 
"30" /im feature in 36 out of 48 C-stars. 

We present in Fig. [I] the IRAS two-colour diagram 
for the sources in our sample following van der Veen & 



Habing (1988 ). There are four sources in our sample with- 



out an entry in the IRAS point source catalogue. For 
these sources we have used ISO/ SWS and LWS obser- 
vations at 12, 25, 60 and 100 /xm to calculate the IRAS 
colours. For IRAS Z02229, no measurements at 60 and 
100 /im are available. In Fig. [I], the warmest sources are 
located in the lower left corner. These are the optically 
visible carbon stars with a low present-day mass-loss rate 
(M ~ 10~ 8 — 10 -7 Mq). With increasing mass loss the 
stars become redder and move up and to the right. After 
the AGB, when the mass loss has terminated, the dust 
moves away from the star and cools; i.e., these sources 
move further to the top-right corner of the diagram. The 
C-stars located above the main group of C-stars have a 
clear 60 /im excess. This is evidence for an additional cool 
dust component. Some of these sources are known to have 



an extended or detached dust shell around them (Young 



Hony et al.: The carrier of the "30" /im feature 



5 



et al. 1993). The empty region between the C-stars and the 
post-AGBs is physical. When the mass loss stops the star 
quickly loses its warmest dust and within a short time span 
(< 1000 yr) the star moves to the right in the two-colour 
diagram. Notice how the sources without a detected "30" 
/im feature cluster on the left of the diagram, i.e., among 
the warmest C-stars. 



2.1. Data reduction 

The SWS data were processed using SWS interactive anal- 
ysis product; IA 3 (see de Graauw ct al. 1996) using calibra- 
tion files and procedures equivalent with pipeline version 
10.1. Further data processing consisted of extensive bad 
data removal primarily to remove the effects of cosmic ray 
hits and rebinning on a fixed resolution wavelength grid. 
If a source has been observed multiple times and these ob- 
servations are of similar quality and of comparable flux- 
level these data are co-added after the pipeline reduction. 
These sources are indicated in Table |l| with a dagger (t). 
Since the features we discuss here are fully resolved in 
all observing modes, we combine the data obtained in all 
different modes to maximise the S/N. Although the wave- 
length coverage of the SWS instrument is well suited to 
study the profile of the "30" /im feature, there are some 
important instrumental effects which hamper the unbi- 
ased extraction of the emission profiles. We discuss these 
below. 



2.1.1. Splicing 

One complete SWS AOT01 spectrum is obtained in 12 
different subbands. These subbands are observed through 
3 different rectangular apertures which range in size from 
14"x20" at the shortest wavelengths to 20"x33" at the 
longest wavelengths. All these data are independently flux 
calibrated and need to be combined to form one continu- 
ous spectrum for one source. We apply scaling factors to 
combine the different subbands to obtain the continuous 
spectra. The C-stars and post-AGB objects we present in 
this study all have a small angular extent even compared 
to the smallest aperture used. Therefore we don't expect 
large jumps to be present due to the differences between 
the apertures used. The angular extent of some PNe can 
be large compared to the sizes of the apertures. If there is 
a clear indication of flux jumps due to aperture changes 
we have not included the source in our sample. 

2.1.2. Leakage 

At wavelengths longer than between 26 and 27.5 /im the 
data of SWS subband 3D are affected by leakage adding 
flux from the 13 /im region. The sources used to derive 
the instrumental response function are all stellar sources 
without circumstellar material. These calibration sources 
are all very blue and emit much more flux at 13 /im relative 
to 26 /im than the cool, red sources we present in this 
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Fig. 2. Examples of the splicing of the SWS band 3D 
(19.5-27.5 /im) and 4 (28.9-45.2 /im) data. We show the 
data before (grey line) and after splicing (black line). All 
data have been scaled to form a continuous spectrum. As 
can be seen; after splicing, the slope of band 3D and band 
4 match. We do not show the band 3E data. The sharp 
rise at 27 /im in R Scl and IRC+40 450 is an instrumental 
artifact (see text for details). 



study. Therefore these calibrators are more affected by 
the leakage than our sources. The instrumental response 
function derived in this way has been implicitly corrected 
for leakage for the blue sources. This resulted in fluxes 
in red sources to be systematically underestimated. More 
recent calibrations (> OLP 10.0), have been corrected for 
this effect. With the improved calibrations, the resulting 
slopes of the spectra beyond 26 /im have been checked and 
are in general agreement with the slope of subband 4. 



2.1.3. The 27.0-27.5 and 27.5-29.0 fim region 

At wavelengths longer than 27.0 /im the data of subband 
3D show a sharp increase which is found throughout the 
complete database of ISO/SWS observations independent 
of source type. The data of subband 3E (27.5-29.0 /im) 
are generally unreliable both in shape and absolute flux 
level. These combined instrumental effects make it inher- 
ently difficult to interpret the 27-29 /im spectra. Any sub- 
structure detected solely in this region alone should be 
distrusted. 
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The instrumental effects between 27 and 29 /im and the 
fact that each of the subbands is independently flux cali- 
brated make it necessary to devise a strategy for splicing 
the band 3D, 3E and 4 data. There is unfortunately no ob- 
jective way to choose this strategy. We choose to assume 
minimal spectral structure between the end of subband 
3D and the beginning of band 4, i.e. to splice the subband 
3D— 4 data in such a way that the matching slopes of 3D 
and 4 also match in flux level. Some examples are shown in 
Fig. ||. The observed discontinuities between subbands are 
relatively small (< 20 per cent) and can be understood as 
the result of absolute flux calibration uncertainties alone. 



2.2. Full spectra 

The resultant spectra for the sources that exhibit a "30" 
/im feature are shown in Fig. |[ The SWS spectra of this 
large group of objects show a spectacular range in colour 
temperature, molecular absorption bands and solid state 
features. The C-stars have molecular absorption bands of 
C 2 H 2 at 3.05, 7-8 and 14 fim, of HCN at 7 and 14 (im, 
CO at 4.7 /im and C3 at 4.8—6 /im. The sharp absorp- 
tion band at 14 is due to C2H2 and HCN. There is an 
emission feature due to solid SiC at 11.4 /im. In the red- 
dest C-stars, we find the SiC in absorption. We also find 
evidence for a weak depression in the 14—22 /im range 
in the reddest objects. This depression could be due to 
aliphatic chain molecules like those found in RAFGL 618 
( |Cernicharo et al. 200l| ). 

The post-AGBs and PNe exhibit many, sometimes 
broad solid state emission features. In many sources we 
find emission due to polycyclic aromatic hydrocarbons in 
the 3-15 /un range. There is a broad plateau feature from 
10—15 /im which may be due to hydrogenated amorphous 
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RU Vir 



carbon (Guillois et al. 1996; Kwok et al. 2001). Many post- 
AGBs and two PNe in the sample have a feature peaking 
at 20.1 /im, called the "21" /im feature in the literature. 
Rece ntly the carrier of this fe ature has been identified with 
TiC (|von Hcldcn et al. 200C| ). The feature at 23 /im found 



in IRAS 18 240 and PN K3-17 is likely due to FeS ( |Hony 



et al. 2002 ). These absorption and emission features have 
to be taken into account when determining the profile of 
the "30" /im feature or the shape of the underlying con- 
tinuum. 

Focusing on the "30" /im feature we can see variations 
in the strength and shape of the band. The most marked 
difference is however a shift in the peak position going 
from 26 /im in some of the AGB stars to 38 /im in the 
PNe. The dashed line in Fig. || and ^ indicates A=26 /im. 
There are systematic changes in the appearance of the 
"30" /im feature from the C-stars to the PNe. The feature 
in the C-stars almost exclusively peaks at 26 /im. There 
are some exceptions like R Scl. In the post- AGB sample, 
the feature is broader and in some sources the feature 
peaks long ward of 26 /im. In the PNe sample, there are no 
sources that peak at 26 /im. However, the appearance of a 
broad feature like the "30" /im feature is sensitive to the 
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Fig. 5. Examples of the fitted continuum. We show the 
spectra (black line), the selected continuum points (dia- 
monds) and the fitted modified blackbody (grey line). 



shape of the underlying dust continuum, especially since 
we have a sample with such a wide range of continuum 
colour temperatures. 




Fig. 3. Overview of the spectra of carbon stars exhibiting the "30" /xm feature. The spectra are ordered according to 
continuum temperature from high to low temperature, bottom to top, left to right. The dashed line marks A=26 fim. 
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Fig. 4. Overview of the spectra of post-AGB objects (left panel) and PNe (right panel) exhibiting the "30" fim feature. 
The spectra are ordered according to continuum temperature from high to low temperature, bottom to top. The dashed 
line marks X—26 fim. The spectrum of RAFGL 618 although warmer than NGC 3918 is shown at the top of the PNe 
for clarity. 
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Table 2. Measured properties. T cont , p are the parameters of the modified blackbody function fitted to the continuum. 
A Ci 3o and P/C are the feature centroid position and peak over continuum ratio. TMgS is the derived temperature of 
the MgS grains. 



cont. "30" /im feature 

Object Tcont p A c , 30 fwhm flux P/C T Mg s 

[K] M [/mi] [W/m 2 ] [K] 



cont. "30" /mi feature 

Object Tcont p A c , 30 fwhm flux P/C T Mg s 

[K] M [ M m] [W/m 2 ] [K] 



NGC 40 
IRAS 00210 
IRAS 01005 
HV Cas 
RAFGL 190 
R Scl 

IRAS Z02229 
RAFGL 341 
IRC+50 096 
IRAS 03313 
U Cam 
RAFGL 618 
W Ori 
IC 418 
V636 Mon 
RAFGL 940 
IRAS 06582 
HD 56126 
CW Leo 
NGC 3918 
RU Vir 
IRAS 13416 
II Lup 
V Crb 
K2-16 

IRAS 16594 
NGC 6369 
IRC+20 326 
CD-49 11554 
HB 5 

RAFGL 5416 
IRC+40 540 



150 
285 0.5 
130 1 
1040 0.2 
275 
2605 -0.2 
235 
380 
855 -0.2 
325 
1775 
235 -1 
2450 
120 1 
1215 
810 
315 
170 
535 
90 1 
1045 
115 1 
625 
1430 
155 0.5 
140 1 
100 1 
770 -0.7 
140 1 
120 
290 
485 



33.6 
28.4 
30.0 
33.5 
30.9 
33.2 
29.1 
29.8 
28.8 
28.6 
31.9 
38.0 
31.3 
30.8 
29.8 
28.2 
29.5 
30.0 
28.6 
33.3 
30.4 
31.6 
29.5 
30.4 
34.4 
29.8 
34.6 
29.1 
30.2 
35.5 
30.4 
28.6 



10.1 

10.7 

11.1 

10.6 

13.0 

13.9 

10.1 

9.4 

9.2 

7.8 

11.8 

10.9 

8.4 

11.3 

10.1 

10.2 

10.3 

12.0 

8.8 

8.5 

10.1 

15.8 

10.1 

10.1 

12.0 

12.1 

10.1 

10.2 

14.0 

11.5 

12.5 

9.1 



5.9e-13 
6.4e-13 
6.6e-13 
1.5e-13 
1.6e-12 
l.le-12 
8.3e-12 
9.4e-13 
1.9e-12 
5.4e-13 
3.9e-13 
5.4e-12 
3.1e-13 
5.5e-12 
1.7e-13 
3.5e-13 
l.le-12 
2.9e-12 
2.7e-10 
7.1e-13 
5.3e-13 
2.8e-12 
3.9e-12 
1.8e-13 
3.4e-13 
9.9e-12 
9.5e-13 
7.4e-12 
4.7e-12 
1.0e-12 
2.2e-12 
8.9e-12 



0.7 
0.8 
1.5 
0.3 
0.3 
1.1 
1.7 
0.4 
0.3 
0.4 
0.6 
0.2 
0.4 
0.9 
0.2 
0.5 
0.4 
0.8 
0.6 
1.0 
0.6 
0.4 
0.3 
0.3 
0.3 
0.9 
1.1 
0.5 
0.7 
0.4 
0.5 
0.6 



110 

300 

220 

100: 

180 

90 

300 

250 

500 

300 

150 

40 a 

150 

180 

250: 

500 

300 

150 

400 

120 

180 

200 a 

400 

150: 

80 

250 

90 

300 

200 a 

70 

220 

400 



T Dra 

RAFGL 2155 
IRAS 18240 
IRC+00 365 
RAFGL 2256 
K3-17 

IRC+10 401 
IRAS 19068 
NGC 6790 
RAFGL 2392 
NGC 6826 
IRAS 19454 
HD 187885 
RAFGL 2477 
IRAS 19584 
IRAS 20000 
VCyg 
NGC 7027 
S Cep 

RAFGL 2688 
RAFGL 2699 
IC 5117 
RAFGL 5625 
IRAS 21489 
SAO 34504 
IRAS 22303 
IRAS 22574 
RAFGL 3068 
RAFGL 3099 
IRAS 23304 
IRAS 23321 



1210 
460 
160 1 
910 



0.3 



390 

100 1 

765 

1165 -0.7 

290 

890 

150 

140 1 

175 

290 

580 

210 

1110 

125 1 

1340 0.1 

200 -1 

540 

130 1 

300 

415 

210 

345 

160 

290 

470 

115 1 

175 



30.2 
28.8 
32.8 
28.6 
29.5 
34.1 
30.0 
28.5 
29.8 
27.7 
32.7 
36.3 
29.6 
30.7 
28.1 
29.4 
30.5 
32.8 
31.2 
31.1 
29.0 
31.2 
30.3 
29.3 
29.1 
30.3 
31.2 
32.4 
29.5 
30.1 
34.5 



10.1 

8.2 

13.1 

11.7 

12.0 

11.5 

10.0 

10.1 

15.6 

8.6 

10.5 

13.1 

10.8 

12.5 

7.5 

12.1 

11.5 

11.0 

9.4 

10.4 

11.4 

9.7 

11.8 

9.7 

10.3 

10.5 

13.6 

14.7 

10.9 

13.4 

13.3 



4.8e-13 
5.7e-12 
1.0e-12 
1.9e-12 
1.9e-12 
1.0e-12 
2.0e-12 
2.0e-13 
9.8e-13 
3.4e-13 
l.le-12 
6.4e-13 
5.2e-12 
2.3e-12 
8.5e-13 
2.5e-12 
1.3e-12 
1.7e-ll 
4.4e-13 
5.9e-ll 
5.9e-13 
7.3e-13 
4.4e-12 
l.le-12 
1.3e-ll 
1.0e-12 
5.9e-13 
8.4e-12 
2.6e-12 
2.3e-12 
6.6e-13 



0.4 
0.6 
1.0 
0.4 
1.0 
0.9 
0.3 
0.4 
1.4 
0.5 
2.0 
0.3 
1.0 
0.6 
1.5 
1.5 
0.3 
0.4 
0.2 
0.3 
0.7 
0.6 
0.4 
0.6 
2.0 
0.7 
0.4 
0.4 
0.7 
1.1 
0.3 
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"Temperature determination uncertain due to optically thick MgS emission. 



3. Continuum 

In order to extract the profiles of the "30" /im features 
and compare them from source to source, we model the 
underlying continuum due to the emission of other circum- 
stellar (CS) dust components. First, we present the way 
we construct these continua and in Sect. |], we discuss the 
resulting profiles. 

To model the underlying continuum we use a simpli- 
fied approach. We represent the continuum with a single 
temperature modified blackbody, 

F(\) = A x B(X,T) x \-p, (1) 



where A is the wavelength, F(X) is the flux density of the 
continuum, B(X, T) is the Planck function of temperature 
T, p is the dust emissivity index and A is a scaling factor. 

We have chosen this approach to estimate the contin- 
uum over doing a radiative transfer calculation for reasons 
of simplicity. The bulk of the CS dust around these sources 
consists of some form of amorphous carbon grains that 
do not exhibit sharp emission features in the wavelength 
range of interest. Therefore, a radiative transfer calcula- 
tion will not yield extra insight into the shape or strength 
of the continuum while introducing many more modelling 
parameters. This method has the advantage that we can 
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compare the feature in such a diverse group of sources in a 
consistent way. Of course Eq. (p]) does not directly allow us 
to incorporate important effects such as optical depth or 
temperature gradients. However varying the p-parameter 
can mimic these effects to some extent. 

The p-parameter reflects the efficiency with which the 
dust grains can emit at wavelengths larger than the grain 
size. Reasonable values of p in the region of interest are be- 
tween 1 and 2. Crystalline materials have this value close 
to 2 and amorphous materials have a p- value between 1 
and 2, while layered materials have an emissivity index 
close to 1. A temperature gradient in the dust shell will 
result in a broader spectral energy distribution (SED). 
This is mimicked by a lower value of p. Likewise an opti- 
cally thick dust shell will result in a broader SED, which 
again can be reproduced by reducing the value of p. 

We use a x 2 fitting procedure to determine the values 
of T and p fitted to selected continuum points in the ranges 
2—22 /tm. If available we also use the LWS spectra to 
verify the continuum at the long wavelength end of the 
"30" /im feature. The 50—100 /mi continuum gives an even 
stronger constraint on the value of p. For most cases the 
resultant continuum runs through the 45 /im region of 
the SWS spectrum. A remarkable exception to this is the 
spectrum of RAFGL 3068. The 2-24 /im spectrum is well 
fitted with a single 290 K Planck function. However we 
find a large excess of this continuum at 45 /xm and the 
available LWS spectrum is not well represented in level or 
slope. Possibly this is due to the optically thick dust shell 
or a biaxial dust/temperature distribution. 

The values for T and p are listed in Table ^. One re- 
markable fact is that the C-stars are well fitted by a sin- 
gle temperature Planck function over the complete wave- 
length range of SWS. The IR SEDs of the post-AGBs and 
PNe are in general less broad and many sources are bet- 
ter fitted with a p- value of 1. We stress however that the 
derived p values cannot be used to constrain the crystal 
structure or the average size of the dust grains in view of 
the aforementioned effects of temperature gradients and 
optical depth. 

4. Profiles 

Using the continua defined in this way, we measure the 
following properties of the "30" /im feature: the centroid 
wavelength (A Cj 3o), i.e., the wavelength where the inte- 
grated flux in the feature at either side is equal; the full 
width at half maximum (fwhm); the flux in the "30" /tm 
feature and the height of the peak of the "30" /im feature 
after continuum subtraction over the continuum ratio, i.e., 
the peak to continuum ratio (P/C). The measured values 
are listed in Table ||. We list upper limits for P/C and 
the flux for the sources without detection. In Fig. ^, we 
show the relation between the [25]- [60] colour and A Cj 3o. 
There is a clear reasonably smooth trend for the feature 
to move to longer wavelengths with redder IRAS colours. 
This indicates that the temperature of the dust is an im- 
portant parameter in determining the profile of the "30" 
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Fig. 6. Centroid of the "30" /im feature with respect to the 
[25]- [60] colour. The symbols are like in Fig. |l|. There is a 
clear trend for the centroid position of the "30" /im feature 
to move to longer wavelengths the redder the object is. 
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Fig. 7. Emissivities for different sources as deduced from 
the ISO spectra. There are large differences in the profile of 
the excess emission. Notice the shift in peak position and 
change in width of the circumstellar "30" /im feature going 
from the C-stars to post-AGBs to PNe (top to bottom). 



/tm feature since [25]- [60] is a direct measure of the dust 
temperature provided that the dust composition in the 
different sources is similar. 

We first try to remove the effect of temperature by di- 
viding by the continuum; a method that is commonly ap- 
plied. Using the modelled continua as described in Sect. || 
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we convert the observed features to relative excess emis- 
sion by dividing by the continuum and subtracting 1. 

If the feature emission is optically thin and the temper- 
ature of the carriers of the feature is equal to the contin- 
uum temperature the derived excess emissions are propor- 
tional to the absorptivity (n a bs) of the carrier and if the 
carrier is the same in these sources then the derived band 
shape should be the same for all sources. However, we find 
large variations in the derived profiles. In Fig. ^, we show 
some examples of the derived profiles. Most notable are 
variations in peak position and the appearance around 26 
/im. Such changes, albeit within a smaller range of fea- 



ture peak positions have led other authors (Volk et al 



2000 , 1 002 ) to conclude that the "30" /jm feature is com- 



posed of two features and the observed variations are due 
to varying relative contributions of these two components. 
One key question is: "What possible causes could there be 
for the observed large variations in band shape?" . We dis- 
cuss three possibilities below. First, optical depth effects. 
Second, temperature effects. Finally, we discuss multiple 
band carriers. 

The optically thin assumption most likely holds be- 
cause the optical depth in the circumstellar shell strongly 
decreases towards longer wavelengths. Note, in this re- 
spect that the "30" /im feature is never found in absorp- 
tion (however, see also Sect. 6.2). Hence, optical depth 
effects are not responsible for the observed profile varia- 
tions. 

Whether the temperature of the amorphous carbon 
grains (defining the shape of the continuum) and the tem- 
perature of the "30" /xm carrier are equal is very uncer- 
tain. The temperature of a dust grain in a circumstellar 
envelope is determined by the distance to the star, the 
absorption properties in the wavelength range where the 
star or the dust shell emits light and the grain size. In 
case the temperature of the grains species responsible for 
the continuum and the "30" /im emission feature are not 
the same, the resulting excess profiles will also not be the 
same from source to source even if the carrier of the band 
is the same. The differences will be very pronounced when 
the emission feature is broad. In this case systematic dif- 
ference between sources are bound to occur in league with 
the strongly changing continuum temperature. Thus, the 
temperature of the carrier of the "30" /im feature is an 
important parameter that determines the profile of the 
emission. 

There may be multiple carriers involved as discussed 
before. In this case the feature near 26 /im dominates in 
the warmest objects while the cooler objects are more and 
more dominated by emission towards 35 /xm. However, this 
scenario has its difficulties since it would require changes 
in the composition of the dust in the relatively dispersed 
and cold nebular surroundings of a post-AGB object or 
even during the PN phase. Such chemical changes can only 
occur extremely slowly, if at all. 

Lastly, variations in grain shape or variations in shape 
distribution can influence the emission profiles. The opti- 
cal properties of materials with a high value of the refrac- 
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Fig. 8. The effect of grain shape. We show the absorp- 
tivity of MgS as a function of grain shape. The numbers 
between brackets refer to the axis ratios of the elliptical 
grains. For comparison we show the derived emissivities 
of RAFGL 190 and NGC 7027. 

tive index are sensitive to the grain shape. Variations in 
the shape distribution will lead to variations in the pro- 
files. 

In our analysis, we will focus on explaining the profile 
variations with temperature variations and the effects of 
variations in the shape distribution of the emitting dust 
grains. 

5. MgS 

Since, it has been demonstrated for a few sources that 



MgS is a viable candidate (Begemann et al. 1994; Jiang 
|et al. 1999| ; |Szczerba et al. 1999| ), we first test MgS as a 



possible candidate for the carrier of the feature. With the 
large sample of good quality spectra in this study we are 
able to test this possible identification systematically in a 
large population of evolved objects. 

As explained above we cannot derive a priori informa- 
tion on the temperature of the "30" /im carrier from the 
observations. Our knowledge is further limited by the fact 
that even for some of the candidate materials like MgS 
or FeS the optical properties are measured only in a lim- 
ited wavelength range. We lack measurements in the UV, 
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Fig. 9. The effect of grain temperature on the MgS emis- 
sion feature. We fold the K a b s of MgS in a continuous dis- 
tribution of ellipsoids (CDE) 

shape distribution with a Planck function of different temper- 
atures. The shape and position of the feature are modified sub- 
stantially. 



shape distribution. The width of the feature calculated us- 
ing CDE matches that of the observed "30" /im feature 
well (e.g. RAFGL 190 in Fig. §, see also Fig. CO). 



5.3. Temperature 

To estimate the MgS temperature (TM g s) we use the con- 
tinuum subtracted spectra with the continuum as derived 
in Sect. |. The emission from MgS grains is calculated us- 
ing the K a b s folded with a Planck function with the tem- 
perature of the grain. Due to the smooth and broad shape 
of the resonance the profile of the emission is very sensitive 
to TMgs- In particular the peak position changes strongly 
with TMgs- This allows us to estimate TM g s from the con- 
tinuum subtracted profiles. This method is most sensitive 
for TMgs < 300 K. Above 300 K, further changes in the 
profile are more subtle since the major part of the feature 
falls in the Rayleigh- Jeans domain of the Planck function. 

We use this temperature estimate and the observed 
band strength in the continuum subtracted spectra to syn- 
thesise a MgS feature in order to compare with the as- 
tronomical spectra. In conclusion, we adopt MgS with a 
CDE shape distribution and allow both the strength and 
the temperature of the MgS grains to vary with respect 
to the underlying continuum. 



optical and near-IR range, which may well dominate the 
dust heating. We have decided to test the MgS identifica- 
tion, leaving the grain temperature as a free parameter. 
We adopt the method we describe below. 

5.1. Material 



We use the optical constants as published by |Bcgemann 



et al. (1994 ). Of the materials they measured, Mgo.9Feo.1S 
is closest to pure MgS. The real and imaginary part of the 
refractive index (n and k values) are given from 10-500 
/im. 

5.2. Shapes 

From the n and k values we can calculate the absorption 
cross-sections for various grain shapes and shape distribu- 
tions in the Rayleigh limit following Bohrcn fc Huffman! 
(1983, chap. 5, 9 and 12). The absorption cross-section of 



MgS around 30 /im is very sensitive to the grain shape. In 
Fig. ||, we show the results of using different grain shapes 
on the absorption properties of MgS. We use a continuous 
distribution of ellipsoids (CDE) for the grain shapes. This 
shape distribution was used by Begemann et al. (1994) 
and gave a good fit to the "30" fim feature observed in 
CW Leo. The same shape distribution was further used 
by |Jiang et al. (1999| ) and [Bzczerba et al. (1999| ). They 
found reasonable fits for the two sources they study. As 
can be clearly seen in Fig. || when comparing the spheres 
with the CDE calculations, the feature broadens and the 
peak position shifts to longer wavelengths using the CDE 



6. Model results 

We show in Fig. [A.l| the observed spectra, the composite 
of the continuum and the synthetic MgS feature and the 
residuals after subtracting the MgS feature. The fits are 
very satisfactory in 50 out of 63 cases. In ~25 sources the 
synthetic spectra obtained with this very simple model are 
able to explain the detailed profile of the "30" /im feature 
very well. The onset and range of the feature and even the 
slight depression between 26—30 /im are reproduced by 
the model. We show a zoomed view of the 30 /im region 
of a few sources that are very well fitted by this simple 
model in Fig. ml. Notice the different apparent shapes of 



the fe ature that the model is able to explain. Volk et al 
(2002|) discuss the "30" /im feature in IRAS 23304 and 
find that they need 2 separate unidentified components 
in order to understand the shape of the feature. Fig. [l(] 
illustrates that this is not necessarily required. 

Examining the complete sample of observed "30" /im 
features and the synthetic spectra, we find there are some 
systematic deviations. In the sample of C-stars and post- 
AGB objects there are numerous examples where the ma- 
jor part of the "30" /im feature is explained well by our 
model, but the observed spectra show excess emission in 
the 26 /xm region. The excess is not accounted for using 
our CDE fits. The most extreme case is IRAS 19584 but 
several sources exhibit the same behaviour. In Sect. 
we discuss the origin of this discrepancy. 

In some cases the synthetic spectra over-predict the 
flux at the longest wavelengths. This can be due to the 
very simplistic method we have used to estimate the con- 
tinuum level. As the dust optical depth decreases with in- 
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Fig. 10. Examples of spectra that are very well fitted with 
MgS in a single temperature, CDE shape distribution. The 
black lines represent the data and the grey line the model. 
The different sources have been offset for clarity. The ex- 



cess ar ound 23 fim in IRAS 1 8240 is due to FeS ( |Kcllcr 
et al. |002| ; [Hony et al. 2002j ). Notice how the model is 
able to explain the profile of the "30" /zm feature found 
in the full range of objects in our sample. 



creasing wavelengths, the continuum level estimated from 
shorter wavelengths might over-predict the true contin- 
uum level. Note, however, that the discrepancies between 
the modelled and the observed spectra in the 26 /im region 
and the 40 /im region cannot be considered completely in- 
dependently. If we were to weaken the strength of the MgS 
feature this would yield a better fit around 45 /im but 
would increase the discrepancy around 26 /im. We also 
note that MgS produces a weak continuum contribution 
at 45 /im (see for example Fig. ^). This continuum con- 
tribution is already taken into account when fitting the 
overall continuum, but it is still present in the calculated 
MgS contribution. Therefore, our model may slightly over- 
predict the fluxes near 45 /im. 

As a class, the spectra of most PNe show another sys- 
tematic difference. The peak position of the "30" /xm fea- 
ture lies in general at longer wavelengths than in the post- 
AGB sources. This is in accordance with the picture of a 
slowly expanding and cooling dusty envelope. We can sim- 
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Fig. 11. Some examples of sources with a 26 /im excess 
compared to the model spectra. We show the data (full 
black line) the model using the CDE calculation (full grey 
line) and the composite model using both spherical MgS 
grains and the CDE calculation (dashed line). We show 
below the emission of spherical MgS grains at a tempera- 
ture of 250 K. The continuum in RAFGL 2155 runs above 
the observed spectrum at the shortest wavelengths, this 



may be due to molecular absorptions (see Sect. 2.2) 



ulate the same shift in peak position using MgS grains. 
However the fits we obtain fail to reproduce the relatively 
narrow width of the observed pro file. We discuss this de- 
viation of the profiles in Sect. |6.3| . 

There are 4 sources in the sample that have a broader 
"30" /im feature than can be fit by our simple model. 
Of these sources, IRAS 13416 and CD-49 11554 show a 
slightly flattened and broadened feature while in the cases 
of RAFGL 618 and RAFGL 2688 the feature is very broad 
with a depression around 30 fim. The latter sources are 
known to have a very large dust column along the line of 
sight. Most likely the feature shape is due to opti cal d epth 
effects. We discuss these sources further in Sect. 6.2. 



Despite these systematic deviations it is clear that our 
simple model is able to explain the profile of the "30" /im 
feature in good detail in a very wide range of objects. We 
conclude that the carrier of the "30" /im feature in the C- 
stars and post-AGB objects is solidly identified with MgS 
and that the variations in peak position reflect differences 
in grain temperature. 

6.1. 26 fim excess 

As discussed above we find about 25 sources that are very 
well fitted by our simple model. Using similar parameters 
we also find about 25 sources which show an excess near 26 
/im. Evidently, there is a contribution from an additional 
dust component in the latter sources. In considering this 
additional dust component we find that MgS itself is the 
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Fig. 12. The effects of optical depth of the profile of the 
MgS emission. We show the absorbed MgS emission fol- 
lowing Eq. (|2|) and (||) (grey dashed and solid lines), the 
profile of the "30" /im feature of RAFGL 2688 (solid black 
line) and the optically thin MgS emission (dashed black 
line) . 



best candidate. The wavelength region where this excess 
occurs (26 /im) is also the wavelength where the generic 
main resonance of MgS occurs. Spherical MgS grains ex- 
hibit a resonance at 26 /im. In our model, the MgS pro- 
file is broader and peaks at a longer wavelength because 
we use the calculated absorption cross-sections in a CDE 
shape distribution. There is no a priori reason why the 
shape distribution should be close to this distribution. In 
this respect, it is rather surprising that our simple model 
works so well for so many sources. We can simulate a dif- 
ferent shape distribution by adding the contributions of 
spherical grains to the CDE profile. This approach is per- 
mitted since we assume optically thin emission and thus 
the contribution of different components add linearly. In 
Fig. O, we show the result of such a composite model 
for one C-star and for one post-AGB object. The relative 
amounts of spherical grains added to the model is ~35 per 
cent. We keep the temperature of the spherical and the 
CDE grains the same as in the initial model. As can be 
seen the spherical grains contribute at the position where 
our initial model fails. We conclude that variations in the 
distribution over grain shapes can explain the variations 
we observe in the profile of the "30" /<m feature. 

6.2. Optically thick shells 

The feature found in CD-49 11554, IRAS 13416, 
RAFGL 618 and RAFGL 2688 is different from the others 
in the sense that it is broader and flatter. This is espe- 
cially true for RAFGL 618 and RAFGL 2688 where the 
profile even shows a central depression. We investigate the 
effect of optical depth on the emission profile. We model 



7i(A) = B(A,T MgS ) x (1 - e-f^ 1 ) and 
MA) - /o(A,T Mg s) x e- p ^ 1 , 



(2) 
(3) 



where TM g s is the temperature of the MgS, B(X, T) is the 
Planck function of temperature T, p is the mass density of 
MgS, I is the column length and Iq(X,T) is MgS emission 
of temperature T. Ii(X) is the limiting case of a column of 
MgS with a single temperature of TivigS- MA) represents 
the case of foreground absorption only; i.e., MgS emission 
of temperature TM g s obscured by a column of MgS with 
negligible emission. In both cases the resulting profile will 
be broader than the optically thin case. For a long column 
Ji(A) approaches B(X,TM g s) while /2(A) becomes double 
peaked with a depression where k\ peaks. The resulting 
ii(A) profile will never show a central depression. We show 
the effects of incorporating the optical depth in Fig. [l2| 
Indeed, these methods yield a broadened profile closer to 
what is observed. The second method reproduces the cen- 
tral depression found in RAFGL 2688. It is important to 
stress that the curves shown in Fig. [l2] are not the result 
of a proper radiative transfer modelling of the CS shell 
or an attempt to fit the observed spectrum of the source. 
Nevertheless, they are able to explain the general char- 
acteristics of the "30" /im profile in the deviant sources 
well. 



6.3. PNe profiles 

For 9 out of 13 PNe the observed "30" /im feature is much 
narrower than in our model. Of the other four cases we dis- 



cussed the profile of RAFGL 618 in Sect. 6.2. It is impor- 



tant to note that the remaining three cases (NGC 6790, 
IRAS 18240 and K 2-16) are the PNe with the highest 
continuum temperature among the PNe in our sample. 

In Fig. ||, we compare the shape of the "30" /tm profile 
of NGC 7027 with the profiles due to differently shaped 
MgS grains. As can be seen an oblate MgS grain with 
an axes ratio of 10:10:1 exhibits a "30" /tm feature which 
peaks at the right position. At present we don't know of a 
physical reason for a preferred oblate grain shape in PNe, 
and a broader CDE shape distribution in the C-stars and 



post-AGB objects (see also Sect. 3.4). 

The shape of a resonance is also influenced by the pres- 
ence of a coating. MgS is very hygroscopic. Under condi- 
tions where oxygen is available in the gas phase MgS can 



be oxidised and transformed into MgO (Nuth et al. 1985 



Bcgemann et al. 1994). It is possible that the MgS is trans- 



formed as the central star of the PN heats up and the UV 
radiation progressively dissociates the CO molecules yield- 
ing gas phase oxygen. This could lead to MgS grains which 
are coated by a thin layer of MgO. We have modelled 
such grains using the electrostatic approximation follow- 
ing Bohren & Huffman (1983, chap. 5). The result is shown 
in Fig. |13|, curve 6. As can be seen the "30" /tm resonance 
is split into two features due to the MgO coating. The fea- 
ture at the red wavelength is shifted to longer wavelengths 
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Fig. 13. The absorptivity of MgS grains coated with vari- 
ous types of materials. The modelled grains are composed 
of 80 per cent MgS and 20 per cent coating by volume. The 
models have been obtained by adding 5 ellipsoidal core- 
mantle grain of axes ratios (10:1:1), (3:1:1), (1:1:1), (3:3:1) 
and (10:10:1) of equal volume. For comparison we show 
the absorptivity of MgS in a CDE distribution (7) and 
the excess in NGC 7027 (8). We show the absorptivity 
of FeS (1) and amorphous carbon (5) grains with a MgS 
coating, both with core volume of 20 per cent of the grain. 
None of the modelled composite grains is able to explain 
the profile of the emission found in many PNe. 



compared to the pure MgS resonance. However the main 
feature is on the blue side of 25 /im towards the strong 
resonance at 18 /jm in the pure MgO material, in clear 
contrast with the observations. 

We explore other possible coatings on MgS grains to 
test their ability to explain the narrow feature observed in 
the PNe and the lack of emission at 26 /im. We find that of 
the composite grains we tested none give a satisfactory ex- 
planation. Mixtures of MgS and FeS have been discussed 



feature ( 


Bcgemann et al. 1994 


Hcnning 2000 


). Curves 1 and 2 



2001 



of embedding an FeS core in a mantle of MgS and em- 
bedding a MgS core in a mantle of FeS, respectively. The 
latter compares most favourably with the position of the 
feature in the PNe. However, the substructure found in 
the spectrum of the composite grain around 33—37 /im is 
not found in the PNe spectra. 

Szczcrba et al. (1999) examine grains of amorphous 
carbon with a mantle of MgS to compare with the 30 /im 



feature in two post-AGB objects. We show simulated spec- 
tra of such grains and MgS grains coated with amorphous 
carbon in Fig. |l^, curves 3 and 4 respectively. Curve 3 
clearly does not match the observed feature in the PNe. 
As can been seen in curve 4 the MgS grains coated with 
amorphous carbon absorb less at 26 /im than pure MgS 
and are therefore a better spectral match to the "30" /im 
feature of the PNe. However, the feature to continuum ra- 
tio in these grains is about a factor 2.5 lower than in the 
pure MgS grains requiring a factor 2.5 more mass in the 
MgS component in order to explain the observed band 
strength. Note also that such grains will still produce a 
weak feature at 26 while in some PNe spectra we find no 
excess at that wavelength at all. 

Lastly, in curve 5 (Fig. FJ|) we show the effect of water 
ice on the MgS grains. The effects on the optical prop- 
erties of a water ice coating are marginal and the profile 
cannot explain the PNe observations. We conclude that 
of the composite materials we have experimented with 
MgS grains coated with amorphous carbon give the best 
spectral match. However we find no composite grains that 
match satisfactorily. 

We stress that although our model does not reproduce 
the "30" /im profile in the PNe in its width it is safe to 
assume that its carrier is MgS based. These PNe are be- 
lieved to be the evolutionary descendants of the sources 
which exhibit the MgS feature. The shift in peak position 
compared to the post-AGB objects follows naturally from 
an expanding and cooling shell. Also, the feature strength 
for the PNe is similar to these found in the post-AGBs 
further strengthening the physical link between the MgS 
in the C-stars and the post-AGBs on one hand and the 
"30" /im feature in the PNe on the other (see also Sect. [7j). 

7. Correlations 

Using the large database of sources available we can study 
some of the properties of the "30" /im feature statisti- 
cally. We have found it most convenient to characterise the 
sources by the temperature of the fitted continuum. We 
use the wavelength where the derived continuum peaks 
(A roax ,cont) as an indicator of the continuum tempera- 
ture. The derived continuum temperature itself is less well 
suited because of the systematic difference in the power 
law index we find between classes of sources. The sources 
are rather uniformly distributed over A max ,cont as well. 

show in Fig. 14(a) the relation between 



we 



First, 

the A max ,cont and the ratio of the integrated flux in the 
"30" /im feature to the total flux in the SWS spectrum 
(Fjo/Isws)- The C-stars demonstrate a clear increase of 
I30/ISWS with decreasing continuum temperature. The 
post-AGB objects emit systematically a larger fraction, of 
up to 25 per cent, in their "30" /im feature. The PNe emit 
a similar fraction in the "30" /im feature as the post-AGB 
objects although with a larger scatter. Notice that the 
sample contains a number of PNe with warm dust indica- 
tive of young PNe. There are a few sources which do not 
follow the general trend. The C-stars, R Scl, IRAS 19584 
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Fig. 14. "30" /im feature properties versus the peak wave- 
length of the continuum. The symbols are the same as in 
Fig. [l| We show in panel (a) the ratio of the integrated flux 
in the "30" /im feature to the integrated flux in the SWS 
spectrum. In panel (b) , we show the peak over continuum 
values. We also show the average values for the C-stars, 
the post-AGBs and the PNe. The non-detections are not 
taken into account in determining the mean values. The 
centroid wavelength of the "30" /im feature is shown in 
panel (c). 
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Fig. 15. The centroid position of the "30" /im feature ver- 
sus the MgS temperature in the model. The symbols are 
the same as in Fig. |l|. We show a power-law function fitted 
to the data with the dashed line. 



and RAFGL 2256, exhibit an atypically strong "30" /im 
feature. These latter two sources are further typified by 
very weak molecular absorptions near 14 /im (see also 
Fig. ||). These observed anomalies are indicative of devi- 
ating conditions in the outflows of these sources, possibly 
a recently halted period of efficient dust formation. The 
post-AGB object IRAS 19454 has a very weak and cold 
"30" /im feature. RAFG L 61 8 has a weak feature due to 
self- absorption (see Sect. 6.2). 

The increasing strength of the "30" /im feature in the 
AGB stars in not surprising. Since the emission is op- 
tically thin I30 is proportional to the amount of MgS. 
The low values of I30/ISWS for the warmest C-stars re- 
flects the fact that there is little dust around these sources 
and most of the IR radiation comes from the stellar pho- 
tosphere. Cooler C-stars have more dust and thus more 
MgS. The difference between the coolest C-stars and the 
post-AGBs is more surprising. The fact that post-AGBs 
emit a larger fraction in the "30" /im feature is due to 
two effects. First since the dust shell becomes optically 
thin in the visible some fraction of the light is emitted at 
shorter wavelengths. Second, the temperature of the MgS 
decreases less rapidly than the temperature of the other 
dust components (see below). 

It is clear that any dust component which produces 
30 per cent of the IR light has to be abundant. In or- 
der to quantify the (relative) amounts of MgS present in 
the CS shells of these objects will require radiative trans- 
fer modelling which is beyond the scope of this paper. 
We can however in first approximation study the relative 
amounts of MgS compared to the other cold dust com- 
ponents by studying the peak to continuum ratio (P/C). 
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In Fig. pl](b), we show the P/C versus the A maX)Con t- The 
majority of the sources lies within the 0.3—1.0 range in 
P/C. We indicate a few clear outliers. R Scl, IRAS 19584 
and RAFGL 2256 have a very strong "30" /im feature in- 
dicating again that these sources have 'too much' MgS for 
a normal C-star. The PNe NGC 6790 and NGC 6826 have 
an exceptionally strong MgS feature. Note that NGC 6790 
also has a very warm continuum, much like a post-AGB 
source or a very young PN. The strong SiC band at 11 /jm 
is consistent with this. We also show the averages for each 
of the classes of sources. The average P/C for C-stars is 
0.5, for post-AGB objects 1.0 and for the PNe it is 0.9. 
The similar ratios for the post-AGB objects and the PNe 
suggests that the carrier of the "30" /im feature in the 
PNe is indeed directly related to the MgS feature in the 
post-AGBs. Furthermore, the similar ranges found for the 
post-AGB objects and the PNe argues against any process 
which results in a destruction of the MgS grains during the 
PN phase. 

In Fig. [l^, we show the derived MgS temperature ver- 
sus the centroid position of the "30" /xm feature. The 
two are well correlated. For convenience, we have fitted 
a power-law function (without physical meaning) to the 
relation. 



MgS 



5.1 10 i4 (A 



-8.34 



c,3Q) 



(4) 



It is not surprising that TjyigS and A Cj 3o are correlated since 
we have used the feature profile to estimate the TMgS- 
However, using Fig. [l5]or Eq. (^) one can easily derive the 
MgS temperature from a given observation. Also indicated 
in the figure are IRAS 13416 and RAFGL 2688; as can be 
seen they fall outside the correlation. This is due to the 
optical depth effects discussed above. 



Lastly, we show in Fig. 16 the relation we find be- 



tween the temperature of the continuum (T cont ) and the 
temperature of the MgS (Tj^gs)- We find, tracing the evo- 
lution from hot dust sources (C-stars) to the PNe, that the 
MgS temperature decreases correspondingly. Surprisingly, 
we find for warmest sources in the sample, with T cont > 
1000 K, very cold MgS. We propose two explanations for 
this phenomenon. First, it may be due to the absorption 
properties of MgS. If MgS cannot efficiently absorb the 
stellar light in the optical or near-IR part of the spectrum 
the grains remain cold. Alternatively, the MgS grains may 
be located further away from the star (see below). In the 
inset, we show a blow-up of the left side of the figure. 
There is a clear correlation between T CO nt and TMgS- The 
MgS in the post-AGB sources is systematically warmer 
than expected on the basis of the C-stars. Apparently, 
in the process of becoming a post-AGB object, when the 
dust shell becomes detached and moves away from the 
star, the general dust cools down more rapidly than does 
the MgS. Note that the difference in dust temperature be- 
tween the coldest C-stars and the warmest post-AGBs is 
seen most clearly in the [12]- [25] colour or A maXiCon t but is 
also found in the [25]- [60] colour (e.g. Fig. |l|). The different 
behaviour of the continuum and the MgS is nicely illus- 
trated in Fig. |l4|(c) where we show that while the contin- 
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Fig. 16. The derived MgS temperature versus the contin- 
uum temperature. The symbols are the same as in Fig. |[ 
We show in the box in the lower right the continuum tem- 
perature of the sources without a "30" /jm feature de- 
tected. The inset shows a blow up on a logarithmic scale 
up of the sources with a continuum temperature below 
1000 K. 



uum becomes redder the position of the "30" /jm feature 
does not change. This effect is the cause for the disconti- 
nuity between the C-stars and the post-AGB object seen 
in Fig. |l4|(a) and (b). The physical cause for this must be 
due to the fact that MgS is actually (partially) heated by 
mid-IR radiation. The mid-IR radiation is much slower to 
respond to the termination of the AGB. Conversely, this 
implies that mid-IR radiation is a less important heating 
agent for the other dust components present. 

8. Discussion 

8.1. Carrier 

The models using MgS we present do very well in explain- 
ing the profile of the feature found in the C-stars and post- 
AGB objects. The model does not do as well in explaining 
the profile of the feature found in the PNe. However, the 
smooth relation found between the [25]- [60] colour and the 
A c , 30 shows a gradually changing profile with dust temper- 
ature. This suggests one carrier or closely related carriers. 
Also the similar P/C values for the different classes of 
objects is indicative of one carrier. 

8.2. The effect of model simplifications 

It is remarkable that our simple model works so well. We 
have applied various simplifications in order to be able to 
study the whole sample in a unified way. These simplifi- 
cations necessarily lead to differences between the models 
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and the observations. Instead of assuming a single tem- 
perature modified blackbody for the shape of the contin- 
uum one can certainly obtain a more detailed fit by tak- 
ing into account the structure of the CS envelope and the 
opacities of the different molecular and dust constituents. 
Our method works well because the contributions of the 
other known other dust components around 30 /im are 
well behaved, smooth and thus indeed form a continuum. 
Moreover, the shape of our modelled continuum between 
20—50 /im is not a strong function of the emissivity in- 
dex (p) and varying the p-parameter or the temperature 
doesn't yield qualitatively different continuum shapes or 
levels. We already pointed out that in some cases our 
model fails at the longest wavelengths of the SWS obser- 
vations, which may be related to our method to estimate 
the continuum. 



We use a single temperature for the MgS to model its 
contribution. In the very extended envelope of a C-star 
or in the nebulous environments in the post-AGBs and 
PNe, the temperature of any dust component will not be 
constant but decreases as a function of the distance to 
the star. The fact that we still get good results using a 
single temperature is a clear indication that the emission 
is optically thin and that the density of the MgS falls off 
sharply with distance. In an optically thin environment 
we know that Td us t oc VR, where R is the distance to 
the star. If the density distribution drops with distance 
as R" 2 or steeper, the contributions are weighted to the 
highest temperature part of the envelope where Xd us t and 
thus the source function is highest. If the density distribu- 
tion is natter there is relatively more dust far away than 
close by. In this case we will observe MgS with a range 
of temperatures. Also sources that are not optically thin 
will emit a feature broader than our single temperature 
MgS model (see also Sect. |6~2| ). We conclude that in the 
majority of the sources the "30" /im emission is due to 
optically thin emission which is dominated by the highest 
temperature MgS closest to the star. 



Recently |Grishko et al. (2001 ) proposed HAC as the 
carrier of the mid-IR emission features in C-stars, post- 
AGBs and PNe. HAC is a very plausible dust compo- 



nent in those environments (e.g. . 


Duley & Williams 1981; 


Goebcl 1987; 


Borghcsi et al. 198 


1: 


Hcnning & Schnaiter 



HAC are dominated by a A 1 continuum (Bussoletti et al 
198 7j). prishko et al. (2001] ) identified some 13 weak spec- 



tral features in the range from 19—120 /im. The strongest 
of which occur at 21,27 and 57 //m. We consider HAC 
a unlikely candidate for the circumstellar "30" /im emis- 
sion feature because the laboratory features are weak in 
contrast to the astronomical data and because the far- 
IR features of HAC are not observed in the astronomical 
spectra. 



8.3. Shape 

We find that a CDE shape distribution fits many observa- 
tions well. A CDE shape distribution is used to model a 
collection of irregularly shaped grains. Irregularly shaped 
grains can for example result from grain growth by ag- 
glomeration ( Blum et al. 200C ). 

We find evidence for differences in the shape distri- 
butions between sources. We have tested for correlations 
between the strength of the 26 /im excess (due to spher- 
ical MgS grain) and the CDE component and other pa- 
rameters like the mass-loss rate, the P/C, the continuum 
temperature or the feature temperature. We find no clear 
correlations. We do however note that we find little evi- 
dence for the 26 /im excess in the hottest C-stars. Stars in 
our sample with a continuum temperature above 1000 K 
do not exhibit the excess. Below 1000 K we find both 
sources with and without the 26 /im excess. Note also 
that the occurrence rate of the 26 /jm excess in the sam- 
ple of post-AGB objects is high. Because the emission in 
the post-AGB phase may be dominated by the dust clos- 
est to the star and hence lost at the tip of the AGB during 
a phase of heavy mass loss (at a rate much higher than 
during the general AGB phase), this might suggest that 
the shape distribution of the grains changes to become 
more spherical towards the end of the AGB, possibly as a 
function of mass-loss rate. 



8.4. Planetary nebulae 

The profile of the "30" /jm feature we find in most PNe is 
narrower than in our model. The largest discrepancy oc- 
curs around 26 /im. Our model fits are able to explain the 
shift in centroid position observed in the PNe; however, 
the models that explain the band-shift best, over-predict 
the flux at 26 and 40 /im. 

There are several observational properties that are im- 
portant when considering explanations for the observed 
discrepancies. First, there is the similar values of P/C in 
the sample of post-AGBs and PNe (Fig. [ljb). This in- 
dicates that the carriers are related and similar in abun- 
dance. Second, the smooth trend we find in the centroid 
position of the "30" /im feature with the [25]- [60] colour 
(Fig. [|). The PNe profiles follow the general trend. This 
indicates that the main effect for the peak shift also in 
the PNe is due to temperature. We also note that there 
are three non-PNe sources in our sample with cold MgS 
(R Scl, IRAS 19454 and IRAS 23321). These sources are 
well fitted by our model. This indicates that the narrow 
profiles are particular to the PN environment. 

We have explored possible MgS based heterogeneous 
grains and variations in MgS grain shape. The heteroge- 
neous grains we have explored do not compare satisfac- 
torily with the observed profile in the PNe. MgS grains 
coated with a layer of amorphous carbon provide a some- 
what better match in terms of the band shape however 
the contrast of the feature with respect to the continuum 
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is strongly reduced. The reduced band strength is in con- 
tradiction with the observed peak over continuum ratios. 

We find that the emission from plate-like MgS grains 
appe ars similar to the feature found in the PNc. In 
Sect. 3T we consider an extra contribution of spherical 
MgS grain to the "30" /jm feature at 26 /im. The vary- 
ing strength of this contribution demonstrates that there 
are variations in the grains shape distribution between 
sources. The question is: "Is it possible that the narrow 
"30" /im feature in the PNe is carried by flattened MgS 
grains?" This would either imply that the MgS grain shape 
distribution in the PNe is heavily skewed to the plate-like 
grains or that the oblate grains emit more efficiently. 

The origin of these variations in the shape distribution 
are unclear. They may reflect variations in the formation 
and destruction of grains. We do not know of a mechanism 
that drives towards oblate grains shapes during the PN- 
phase or a mechanism that selectively destroys spherical 
and prolate grains. Note, that the peak over continuum 
value in the PNe is similar to those found in the post- 
AGB objects arguing against a destruction of a part of the 
MgS grains. It is important to remark that in the CDE 
distribution only ~20 per cent of the grains is oblate with 
an axis ratio of l:t or more extreme, i.e. with axes ratios of 
1:X:Y, where Y is smaller than ^ and X can vary between 



' X and 1 . Such axes ratios are required in order to shift 
the peak position to the wavelength where the "30" /jm 
feature in the PNe peaks and to suppress the emission at 
26 /im. 

Alternatively, these inferred shape variations may re- 
sult from a shape dependent heating of grains. In the 
C-stars, likely, heating and cooling occurs at IR wave- 
lengths. Because the grains are small compared to these 
wavelengths the temperature will be shape independent. 
However, in the PNe heating occurs through absorption 
of visible and UV radiation, which scales with the cross- 
section. Hence in that case dust temperatures will be 
shape dependent. Indeed flattened grains have on average 
a larger geometric cross-section compared to their volume 
that spherical or prolate grains. 

As noted before, the PNe in the sample with the 
warmest dust continuum are well fitted by our model. 
The model fails in the PNc sample at MgS temperatures 
below ~90 K. Perhaps, the observed discrepancies are 
due to changes in the optical properties of pure MgS at 
low temperatures. No laboratory measurements of MgS 
at these temperatures are published. FeO has the same 
lattice structure as MgS and measurements of FeO at dif- 
ferent te mperatures (10 K, 100 K, 20 K and 300 K) are 
available (rlcnning & Mutschke 1997). We have used FeO 
as an analogue to examine the effect of temperature. The 
trend with lower temperature is for the FeO resonance at 
20 /im to become narrower and stronger with respect to 
the continuum. There is a small shift in the peak position. 
Comparing the 300 K measurement to 10 K the latter 
resonance peaks ~0.1 /im more to the blue for spherical 
grains. Like MgS, FeO is also very sensitive to the shape 



of the grains. Using a CDE shape distribution the reso- 
nance of the 10 K sample lies ~1 /im on the blue side of the 
300 K resonance. If we translate this behaviour to the MgS 
data and the "30" /im feature in the PNe it worsens the 
situation since we would expect a stronger contribution 
at shorter wavelengths exactly where our model already 
over-predicts. As an alternative, one might speculate that 
a change in lattice structure at low temperatures may oc- 
cur resulting in different optical properties. It is interest- 
ing to note that Bcrthold (1964) reports that MgS can 
condense in two lattice structures: cubic and hexagonal. 
This author finds that hexagonal MgS exhibits a narrower 
mid-IR resonance than cubic and amorphous MgS. 

We conclude that, although it is possible to find MgS 
based candidate materials that provide a better spectral 
match to the observed "30" /im feature in the PNe no ex- 
planation for the "30" /im profiles in the PNe is completely 
satisfactory at this time. 

8.5. MgS Temperature 

We have found that the "30" /im feature is very well ex- 
plained with MgS provided that the continuum tempera- 
ture and the MgS can be substantially different from each 
other (see also Fig. Moreover, the temperature dif- 
ference changes between classes of sources. In the C-stars 
we find that MgS is colder than the continuum while in 
the post-AGB sample the MgS grains are warmer. The 
temperature difference and the fact that the continuum 
temperature varies relative to the MgS temperature im- 
plies that the MgS grain are not in thermal contact with 
the amorphous carbons grains that carry the continuum. 
Thus, MgS is condensed in these environments into sepa- 
rate grains. 

The relatively warm hot MgS that we observe in the 
post-AGB sample is certainly a phenomenon that deserves 
further study. As pointed out above, this must be due to 
a difference in the heating properties of MgS relative to 
the other dust constituents. We suggest that due to the 
strong resonance in the mid-IR MgS is (partially) heated 
by IR radiation. This may cause a temperature difference 
as observed. It requires measured optical constants for 
MgS over a much wider wavelength range than now avail- 
able and radiative transfer modelling to understand the 
detailed temperature behaviour of the MgS grains. 

We have shown in Fig. [l(] that the warmest C-stars 
that exhibit the "30" /im feature have cold MgS. We con- 
sider two possible explanations for this behaviour. First, 
MgS can be cold because it is unable to absorb the radia- 
tion from the central source efficiently. This would require 
that the MgS grains are very transparent in the visible 
and near-IR where these warm C-stars emit most of their 
radiation. In this case a cooler star will emit more radia- 
tion at wavelengths where MgS can absorb and as a result 
the MgS grains will be warmer. This explanation is con- 
sistent with our notion that MgS is (partially) heated by 
IR radiation. 
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A second intriguing possibility is that the MgS around 
these sources is located far away from the star. R Scl, 
U Cam and W Ori are known to have a detached dust 



shell (Zuckerman 1993). These shells are formed during 
an earlier phase when the mass-loss rate was higher than 
the current mass- loss rate ( |Willems fc de Jong 1988 ). If 
during this phase the AGB star was already a C-star, MgS 
could have condensed in the outflow and be present in the 
detached shell. We have sketched in Fig. |l| the evolution of 
a star which suffers a brief period of enhanced mass loss. 
First, the star becomes redder during the phase of high 
mass loss. When the shell becomes detached the 12 /mi 
flux quickly drops while the inner edge of the dust shell 
moves away from the star and the warmest dust is rapidly 
lost. The star moves back to the locus of the warmest 
C-stars but with an excess of cold dust. The excess of 
cold dust is observable as a 60 /im excess. The sources 
mentioned above indeed show the 60 /im excess. We note 
that some stars with 60 /im excess do not show the "30" 
/im feature. 

To distinguish between these two possibilities requires 
radiative transfer modelling and a more detailed investiga- 
tion of these particular sources which is beyond the scope 
of this paper. The location and temperature of the MgS 
in the warmest C-stars will be further discussed in (Hony 
et al., in prep.). 

8.6. MgS in the ISM 

The behaviour of the MgS we find in this study has im- 
plications for our understanding of the fate of MgS that 
is produced during the C-star phase. Previously, the exis- 
tence of MgS in the ISM has been excluded because of the 
lack of spectral signature at 26 /im. We have shown that at 
low dust temperatures no emission signature is expected 
at 26 /im. For the dust temperatures in the ISM we would 
expect a smooth feature at wavelengths longer than 35 
/im. In this respect it is interesting to note that the peak 
to continuum (P / C) values in the post- AGB and PNe sam- 
ple are similar. We find no evidence for a rapid destruction 
by the UV radiation field in the PN. This shows that the 
MgS will indeed be injected into the ISM. Of course, in the 
ISM the contrast will be much less because the injected 
material will be diluted with other star dust and general 
ISM material. Nevertheless a search for the "30" /im fea- 
ture in the ISM may be very worthwhile. In this respect it 
is important to note that Chan et al. (1997 ) have discussed 
a broad emission feature they observed along several lines 
of sight towards the galactic centre. The feature they ob- 
served is similar to the "30" /im feature observed in the 
C-rich post-AGB objects. 



9. Summary and concluding remarks 

We present a large sample of 63 ISO/SWS spectra of 
carbon-rich evolved objects that exhibit a "30" /im fea- 
ture. The feature is detected along the complete evolution- 



ary track from low mass loss C-stars and extreme C-stars 
via post AGB objects to PNe. 

We present a simple approach to determine the con- 
tinuum. We study the profile of the "30" /im feature after 
continuum subtraction. We find large systematic varia- 
tions in the appearance of the "30" /im feature. We firmly 
identify the carrier of the feature with MgS in a CDE 
shape distribution. The profile of the observed feature is 
consistent with MgS provided that the temperature of the 
MgS can be different from the bulk of the dust. This ap- 
proach allows us to model, in a unified way, the profiles of 
the "30" /im feature in a wide range of objects even when 
the feature can appear extremely different. 

We find an additional component at 26 /im in ~25 
sources. We argue that this component is due to differ- 
ences in the distribution over shapes of the MgS grains; 
specifically, it requires a distribution more weighted to 
spherical MgS grains. We find no clear correlations of 
this excess with other properties of the sources. The 
self-absorbed "30" /im features of RAFGL 618 and 
RAFGL 2688 reflect the high optical depth in the "30" 
/im feature in these sources. 

We find that the typical profile of the "30" /im fea- 
ture in the PNe is narrower than predicted by our model. 
We consider several possible explanations. We find that 
flattened MgS grains provide a better spectral match to 
the "30" /xm feature in the PNe. However it is presently 
unclear why in the PNe environment the "30" /im feature 
is dominated by flattened grains. 

We find that the temperature of the MgS grains is 
different from the bulk of the dust. Therefore, they can- 
not be in thermal contact with the other dust species but 
must exist as separate grains. In the C-stars MgS is colder 
than the other dust while in the post-AGB objects MgS 
is warmer. Likely, this is because the MgS is efficiently 
heated by mid-IR radiation which is less important in the 
C-stars. 

The behaviour of the temperature of the MgS for the 
hottest C-stars is enigmatic. The hottest sources have very 
cold MgS grains. We propose two explanations for this 
phenomenon. First, the "30" /im feature in these sources 
may be due to a previous mass-loss phase and the MgS 
is thus located far from the source and cold. Second, the 
MgS grains absorb very inefficiently in the optical and 
near-IR and therefore the hottest C-stars do not heat the 
MgS grains well. 

We examine the feature over continuum ratio to study 
the relative proportion of MgS to the other dust compo- 
nents. We find no evidence for rapid destruction of MgS 
grains during the PNe phase and possibly MgS grains may 
survive to be incorporated into the ISM. 

We would like to emphasise the need for further mea- 
surements of the optical properties of astrophysically rel- 
evant materials at all relevant wavelengths. This study 
has put forwards some very interesting questions in rela- 
tion to the absorption properties of MgS in the optical 
and near-IR and also concerning low temperature effects. 



Hony et al.: The carrier of the "30" /j,m feature 



21 



Unfortunately, no published laboratory measurements are 
available to test some of our scenarios and conclusions. 
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Fig. A.l. "30" /jm feature modelled using MgS optical properties. 




Fig. A.l. (continued). 
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1. (continued). 
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Fig. A.l. (continued). 
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1. (continued). 
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Fig. A.l. (continued). 
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Fig. A.l. (continued). 
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Fig. A. 2. C-stars with a weak or without a "30" /im feature. We show the fitted continuum (dashed line). The sources 
are sorted according to decreasing continuum temperature. We refrain from showing a MgS model. In those cases 
where the strength of the "30" /im feature allowed to extract a temperature of the MgS the value is listed in Table |[ 
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Fig. A. 2. (continued). 



